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Abstract. The dynamical response of a helmet streamer to a flux rope escape from the sub-photosphere 
is examined in a physically self-consistent manner within the approximation of axisymmetric three- 
dimensional magnetohydrodynamics (i.e., so-called *2j D’). In contrast to the previous planar ana- 
lyses of Paper 1 (Wu, Guo, and Wang, 1995), the present study shows, with the inclusion of out- 
of-plane components of magnetic and velocity Helds, that the magnetic configuration represents a 
helical flux rope instead of a planar bubble as shown in Paper 1. Because of this more physically- 
realistic configuration, we are able to examine the dynamical evolution of the helical flux rope’s 
interaction with the helmet streamer. This process leads to the formation of two parts of the solar 
maw ejection; (i) the expulsion of the helmet dome due to eruption of this flux rope, and (ii) the 
flux rope’s eruption itself. When this two-part feature propagates out to the interplanetary space, it 
exhibits all the physical characteristics of observed interplanetary magnetic clouds. These numerical 
simulations also show that the dynamical behavior of the streamer-flux rope system has three distinct 
states: (i) quasi-equilibrium, (ii) non-equilibrium, and (iii) eruptive state depending on the energy 
level of the flux rope. 


1. Introduction 

In the first paper of this series (Wu, Guo, and Wang, 1995, hereafter Paper I), we 
presented a two-dimensional planar magnetohydrodynamic simulation to invest- 
igate the dynamical evolution of a coronal helmet streamer due to emerging flux 
in the form of a magnetic bubble, namely, a detached magnetic structure that has 
been introduced into the closed region of the streamer. The dynamical behavior 
of this magnetic structure is determined by the interaction of the three forces, 
(i.e., magnetic force, gravity and thermal pressure gradient) (see Figures 8-9 of 
Paper I). The numerical simulation shows three distinct states according to the 
strength of this detached structure. These three distinct states are: (i) equilibrium, 
(ii) non-equilibrium, and (iii) eruptive states; these states enable us to understand 
the physical processes of destabilization of coronal helmet streamers. Recently, 
Hundhausen (1993) suggested that most observed non-flare associated (loop-like) 
coronal mass ejections (CMEs) are caused by destabilization of coronal helmet 
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streamers. The density features of this and earlier planar MHD simulations did not 
provide resemblance to the observed loop-like CMEs (Sime, MacQueen, and Hund- 
hausen, 1984). This is, we believe, due to the limitation of the two-dimensional 
model in which the B $ component of the coronal magnetic field, B, could not 
be included. The objective of the present study is to improve this MHD simula- 
tion model and, thereby, to enable the results to exhibit typical observed loop-like 
CMEs. We have extended our previous two-dimensional planar MHD model in 
pursuit of this objective to the three-dimensional axisymmetrical case (i.e., so- 
called ^-dimensional model’). In this model, the out-of-plane components of 
the magnetic field and velocity are included; therefore, the magnetic bubble in the 
planar model becomes a flux rope in the present model as shown in Figure 1. 

Recently, a number of authors (Chen and Garten, 1993, Chen 1989; Kumar 
and Rust, 1996) have investigated the dynamics of the flux rope in the sol- 
ar/interplanetary space. Because of the importance of flux-rope dynamics as a 
cause of interplanetary disturbances, such as interplanetary magnetic clouds (Bur- 
laga, Behannon, and Klein, 1987; Lepping, Jones, and Burlaga, 1990), they lead 
to geomagnetic activity if the Earth is intercepted (Wright and McNamara, 1983; 
Tsurutani et al. 1988; Gosling et al., 1991, 1994). However, these studies by the 
above-noted workers are in a form of an analytical approach, thus a complete 
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MHD description in a self-consistent manner could not be included. In addition, 
the background of the undisturbed solar-interplanetary atmosphere is treated in an 
ad-hoc basis. In the present study, we will present a self-consistent MHD model 
in D with a background atmosphere that has a helmet streamer and a flux-rope 
under the helmet dome. Then, we will watch this system evolve dynamically on the 
basis of MHD theory. Section 2 presents the mathematical model and procedures 
to construct the helmet-streamer and flux-rope system. The numerical results are 
included in Section 3. Physical interpretations and concluding remarks are given 
in Section 4. 


2. Mathematical Model and Method of Treatment 


2.1. Mathematical model 


The basic equations appropriate to the present study are in the system of spherical 
coordinates (r, 9 , <p), where r is the radius, 9 the co-latitude, and ip the longitude 
with the assumption of axisymmetry (i.e., d/d(p = 0). Thus, the equations of mass, 
momentum, and energy conservation, together with the induction equation, in the 
ideal MHD approximation are as follows: 
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The parameters are: u r the radial velocity; the meridional velocity; u v the 
azimuthal velocity; B r , Be, and B v the components of magnetic field; 7 the 
polytropic index (7 — 1.05); po the magnetic permeability; G tje gravitational 1/ 
constant; M s the solar mass; p the thermal pressure; p the density; and T the 
temperature. The equation of state, p = 2 nkT, is used with n as number density; 
and k as the Boltzmann constant. 


u e 

the. 


2.2. Method of treatment 


The method of treatment used in this paper is similar to that employed in Paper I 
to which the interested reader is referred for the details. However, the initial state 
of this study is a helical flux rope in 2^-D nonplanar geometry imbedded in the 
closed field region (i.e., helmet dome) of a coronal helmet streamer as shown in 
Figure 2(b) which is in contrast to the streamer-bubble in a 2-D planar configura- 
tion used in Paper I. Therefore, the flux rope is approximated by an axisymmetric 
magnetic toroid solution given by Shaffanov (1960). This solution describes an 
equilibrium magnetic toroid which satisfies the Grad-Shafranov equa^on in the 
approximation of small curvature (see Shafranov, 1960, for details). The expres- 
sions of this solution for plasma pressure (p) and magnetic field (B) in cylindrical 
coordinates are: 
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(a) 




(c) 


Figure 2. (a) The initial magnetic field configuration of the helmet streamer and flux-rope, (b) The 
computed magnetic field configuration of the streamer and flux-rope system after the flux-rope 
emerged into the closed field region of the streamer for Case 1. (c) The schematic description of the 
three-dimensional view of the streamer arcade system with a filament in it. 
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current (see Shafranov, 1960, for details). The radius, pressure, magnetic field and 
currents are normalized, respectively, as follows: 


R! — r/R a , 
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(6) 


with , p s , and B s being solar radius, plasma pressure, and magnetic field strength 

on the solar surface, respectively. For convenience, the prime is dropped in Equa- 
tions (2)-(5) that are the final expressions for an axisymmetric magnetic toroid 
in cylindrical coordinates. In order to introduce these expressions into the present 
simulation, we have to transfer them into spherical coordinates (r,0, <p)as shown 
in Figure 1 and the parameters to define the flux rope are shown in Table I. For 
convenience, constants and J v are replaced by and which have clearer 
physical meaning. We set R — R s in the torus solution, thus, only four parameters 
remain to determine the solution. It is worth pointing out that Equations (3) — (5) 
analytically satisfy the solenoidal condition (V • B = 0). In order to assure the 
numerical accuracy of the solenoidal condition, the divergence-cleaning procedure 
is again implemented in the present calculation (Ramshv^ 1983) as we did in 
Paper I. 

To proceed with this simulation, we follow the same procedure as described in 
Paper I by first placing the flux rope below the streamer as shown in Figure 2(a), 
then allowing the flux rope to move very slowly with respect to the Alfven speed 


/ 0 
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Table I 

Flux rope parameters 


Flux rope radius 


nrf'V_ / 


R Large radius (the center of the torus refers to the center of the Sun 
po Thermal pressure at the surface of the torus 

f3\ plasma (3 at the central axis of the torus 

fa plasma P at the surface of the torus 


Table n 

The three sets of parameters for the descrip- 
tion of the flux-rope 


Case 

a 

P0 

0i 

0i 

1 

0.039 R, 

4.0 p a 

1.0 

0.5 

2 

0.045 R, 

4.0 p. 

1.0 

0.5 

3 

0.1 R, 

0.75 p. 

0.2 

0.1 


into the computational domain. The solution that represents the helmet streamer 
is identical to the one given in Paper I. After the flux-rope moves into the closed 
field region of the streamer, the numerically relaxed solution for Case I which 
represent the streamer and flux-rope system is shown in Figure 2(b). On the other 
two cases, they have similar configuration but are at transient state. These cases can 
be interpreted as a streamer-arcade with a filament in it, as shown in Figure 2(c). 


3. Numerical Results 

The numerical simulation was carried out for three sets of parameters for the flux- 
rope (i.e., toms). These three sets (or cases) of parameters are given in Table II. 

Using these three sets of representative physical parameters, we again obtained 
three numerical solutions with distinct states as found in the two-dimensional 
planar case shown in Paper I. These three distinct states could be classified into 
three categories; (i) quasi-equilibrium, (ii) non-equilibrium, and (iii) eruptive state. 
These solutions are shown in Figure 3. The amount of total energy of the flux-rope 
before entering the streamer in each of these three cases is shown in Table HI. It 
is readily seen from Table m that the eruptive state has the highest total energy 
content which includes the thermal and magnetic energy. This case is followed by 
the non-equilibrium and quasi-equilibrium states, respectively. It is also interesting 
to note from among these three cases, that the magnetic energy increased 75 times 
from the quasi-equilibrium state to the eruptive state, whereas, the thermal energy 
increased only by a factor of 16. The latter thermal energy increase is due to 
the increased mass content in the flux-rope. This comparison indicates that the 
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Figure 3. Radial distance of the loci of the center of the flux-rope vs time for three cases; (1) quasi- 
equilibrium state, (2) non-equilibrium state, and (3) eruptive state. 


Table in 

Energy content of 2j-D magnetic flux-rope (Torus) 


Case 

Thermal energy* 

Magnetic energy* 


(«gs) 

(ergs) 

1 (o = 0.039 R a ) 

3.61 x 10 30 

2.33 x 10 29 

2 (a- 0.045 R , ) 

4.80 x 10 30 

3.10 x 10 29 

3 (o = 0.100 rt 4 ) 

5.71 x 10 31 

1.75 x 10 31 


®To compute the energy content, we have assumed the length ¥■ 

of the flux rope (i.e., third dimension) as 0.1 


magnetic energy increase may be the primary source and responsible physical 
agent for the eruption. This can also be seen from an analysis of the interaction of 
the three forces as shown in Figure 4 and discussed below. 

We select the eruptive case (i.e.. Case 3) for this discussion to illustrate how this 
flux-rope destabilizes the helmet streamer to form the CME and its interplanetary 
consequence (i.e., observed interplanetary magnetic cloud). 


3.1. Initiation of solar mass ejections (SME) 


As we pointed out in Paper I, the acceleration of the magnetic bubble in the 
two-dimensional planar case, which may be considered as a flux-rope without 
twist, is due to the nonlinear interactions of the three major forces (i.e., pressure 
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THREE FORCES VS R CURVE (T - 6 HOURS) THREE FORCES VS R CURVE (T - 6 HOURS) 




Figure 4. Distributions of the pressure gradient force ( — ), magnetic force ( ), and gravitational 

force ( ) normalized to the first term of the momentum equation at the equator (left panels) and 5° 

away from the equator ( right panels) for time 4 and 6 hours after flux-rope emergence, respectively, 
for Case 3 (see Table HI). 


gradient, magnetic, and gravitational forces). Figure 4 shows these three major force 
distributions at the equator and five degrees off the equator in the left and right 
panels, respectively, four hours and six hours after the emergence of the flux-rope 
into the solar atmosphere for Case 3. Figure 4 shows that behind the high density 
loop (represented by the outward most peak of pressure gradient), there is always 
a peak of magnetic force which can be interpreted as the magnetic force pushing 
the loop outward, however, closer to the solar surface, there is a complex interplay 
among the three forces (i.e., magnetic, pressure, and gravitational forces). The loop^ 
itself is clearly driven by the outward Lorentz force (J x B) which is represented 
by the dotted line. This result implies that the propagation and acceleration of the 
flux-rope is magnetically driven. 

Figure 5 shows the evolution of magnetic field lines, velocity vectors, and density 
contours at 2, 3, and 4 hours after the emergence of the flux-rope. Corresponding 
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HELD UNC AMO VELOCITY POD (T- 2 HOURS) 


contours or ocxsity (t- 2 hours) 


nao uc ano velocity nao (r- 3 hours) contours or ochstty (t» j hours) 



7 8 S a 3 2 T 1 2 3 4 3 R 7 7 8 S A 3 2 < T 2 3 4 5 

SOUR RAM SOLAR RAO* 



mo urc amo velocity reuj p- * hours) contours or ochstty (t- a hours) 



SOLAR RAH SOUR RAO* 

Figure 5. The magnetic field lines, velocity vectors, and density contours for Case 3 at t — 2, 3, and 
4 hours after flux-rope emergence. 

plots in a three-dimensional representation are shown in Figure 6 for the magnetic 
field lines and relative polarization brightness, pB. The latter is computed from 
an integration of the density along the line of sight. From these results, we notice 
that there are two parts for this solar mass ejection. The first part is due to upward 
motion of the flux-rope which destabilizes the streamer and, thereby, causes the 
coronal mass in the helmet dome to form the bright loop of the usually observed 
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loop-like CMEs (Figures 5 and 6). These results indicate that the greatest density 
enhancement occurred at the flanks of the loop rather than at the loop apex. This 
result is similar to the numerical 2D MHD model results given by Steinolfson and 
Hundhausen (1988) who used a pressure pulse and volumetric heating to generate 
aCME. The second part of the brightness is caused by the material contained in the 
flux-rope which may be considered to be the erupted filament. It is worthy to note 
that the right side panel of Figure 6 shows the 3-D magnetic field configurations 
of the streamer and flux-rope. At 2 hours, the flux-rope is hidden behind the 
occulting disk, which is shown as a dark circle. This procedure is necessary for 
comparisons with coronagraph observations when we need to use the vignetting 
function for integration of the density along the line of sight in order to present 
the relative polarization brightness. At the times of 4 and 6 hours, we clearly show 
the propagation of the flux-rope to interplanetary space, thereby leading to the 
formation and, presumably, detection of the magnetic cloud. A detailed discussion 
of this feature is given in the next section. 

3.2. Propagation of the flux rope 

We extended our computation to 30 R a in order to examine the interplanetary 
consequences of this streamer and flux-rope system. A grid system of 210 x 62 is 
used for this purpose. Figure 7 shows the magnetic field lines, velocity vectors, and 
density contours at 8 and 12 hours, respectively, after the emergence of the flux- 
rope for Case 3 (i.e., eruptive case). Figure 7 is, in fact, the continuation of Figure 6. 
By examining the evolution of magnetic field lines as shown in Figures 6 and 7, we 
recognize that the flux-rope has ballooned out, pushing the helmet dome laterally 
aside, and developing a magnetohydrodynamic fast shock as shown in Figure 8(a). 
In Figure 8(b) we present the distribution of the MHD fast wave Mach number. 
Each point on the curve in Figure 8(b) represents the position of the wave front and 
its corresponding fast wave Mach number at a specific time. This figure shows that 
the MHD fast shock begins to develop at ~ 2.5 R a for this case. This fast MHD 
shock is caused by the movement of the flux-rope as a piston. It should be noted 
that in Figure 8(b), the fast wave Mach number does not start at 0. This is because 
once we start moving the flux rope into the computational domain, it will generate 
a fast MHD wave. At the very early stage, the real driver, i.e., the flux rope which is 
expanding slowly, is behind the wave and is not distinguishable. Once the flux rope 
begins to accelerate, it will generate another wave that will eventually evolve into a 
fast MHD shock that overtakes the foregoing MHD wave. The radial velocity of the 
center of the flux-rope is plotted for both the eruptive and non-equilibrium cases 
as seen in Figure 9. We observe that the radial velocities reach asymptotic values 
of 420 km s" 1 and 380 km s -1 , respectively, for these two cases. The strengths of 
the azimuthal component of the magnetic field, corresponding to these two values, 
are 7.24 and 6.32 G, respectively. This means that the shear angle of the twist of 


solah711.tex; 15/10/1996; 9:15; v.6; p.ll 






HELMET STRAMER AND FLUX-ROPE INTERACTION, II 


13 


mo unc «o mom nan cr- s hours) 


CONTOURS V BOOTY <T- 6 HOURS) 



Tsnntiias » m it aa a x » t» 10 t» it i i w it a a a 

wuoconwe oaxmct <*•) tRueccniNC c awcr <*») ■ 


mo unc ano voocar mo (r- 12 hours) contours of KMnr {?• 12 hours) 



■ » a it r 1 1 uuaaa »»»i»»it 1 a 11 a a a 


HOJOODRNB OSOMOC (Da) WUOCBRNO OBNCC (*») 

Figure 7. The magnetic field lines, velocity vectors, and density contours for Cases 3 at 8 and 
12 hours, respectively, after the flux-rope emergence. 


the flux rope is the key parameter in the acceleration process for the development 
of shocks. 

Finally we examine the evolution of this streamer flux-rope system in the inter- 
planetary space. We plot, in Figure 10, the density, temperature, radial velocity, and 
magnetic field as a function of time, at 25 R s and a displacement of 15° away from 
the equator for the eruptive case. The magnetic field in Figure 10 is represented 
in solar ecliptic coordinates in order to compare with observations conveniently at 
both the hypothetical spacecraft location as well as at larger distances. If we com- 
pare Figure 10 with the observations of magnetic clouds at Earth (e.g., Burlaga, 
Behannon, and Klein, 1987), the local density and temperature minima behind 
the shock are noted, thereby immediately indicating that they are qualitatively the 
same. Quantitatively, they are certainly different, because the present results are 
at 25 R s . We perform several simple extrapolations using the polytropic law. For 
example, the corresponding temperature minima at 1 AU are: 2.2 x 10 5 K for 
7 = 1.20; 4.4 x 10 4 K for 7 = 1.4; and 0.5 x 10 4 K for 7 = 1.67. From these 
results, we may conclude that when the solar wind expands from the corona to the 
interplanetary space, additional heating is needed (as expected) for both transient 
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DISTANCE (SOLAR RADII) 


Figure 9. The velocity distribution of the flux-rope center for Cases 2 and 3, respectively, corres- 
ponding to the non-equilibrium and eruptive cases. 


as well as quiescent flows (Suess, Wang, and Wu, 1996). From the density plot of 
Figure 10, we identify two peaks; the sharp one is the fast MHD shock which is 
caused by the propagation of the helmet dome, and the fatter one corresponds to 
the flux-rope. These features are related to our suggestion of solar mass ejections 
that consist of two parts. 


4. Physical Interpretations and Concluding Remarks 

The self-consistent planar MHD simulation model for the study of the dynamical 
evolution of a coronal streamer-bubble (Paper I) has now been extended to the 
nonplanar case (i.e., three-dimensional, axisymmetric geometry). This extension 
made it possible for us to study the dynamical evolution of a helical structure 
within a helmet streamer. Physically, this configuration corresponds to a flux- 
rope emerging from subphotosphere into the helmet dome of a coronal streamer. 
Thus, we used this model to examine the outward motion of this flux-rope and 
its interaction with a helmet streamer. The results obtained from this simulation 
demonstrated, clearly, that there are three distinct states; (i) quasi-equilibrium, 
(ii) non-equilibrium, and (iii) eruptive state as shown in Figure 3. These three 
numerical solutions could be interpreted as three different physical states of a flux- 
rope that may be embedded within the solar atmosphere. If the energy content of the 
flux-rope (i.e., B v is small) does not reach a threshold, the energy will remain within 
the flux-rope in a stationary position for a long time; the consequence is formation 
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Figure 10. The density, temperature, radial velocity, and magnetic held variations versus time at 
r — 25 R a and 15° (away) from the equator for Case 3. The magnetic fields arc given in solar ecliptic 
coordinates (Burlaga, 1988) for direct comparison with a hypothetical spacecraft at that location. 

of a stationary filament (i.e., Figure 3, Case 1). If the filament reaches the threshold, 
its short term future will depend on how far the filament energy content exceeds the 
threshold. That is, two characteristics are exhibited; one shows a period of a few 
hours agitation, then outward motion, slowly accelerating over several hours, then 
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reaching an asymptotic speed propagating outward; the other one is the flux-rope’s 
£ immediate emption. These two characteristics are represented by the solutions 
shown in Figure 3 (i.e., Cases 2 and 3). When this flux-rope takes off, it compresses 
the helmet dome and pushes it away. These properties are shown in Figures 5 and 
6. It is clearly recognized that this eruptive flux-rope has pushed away the coronal 
mass and magnetic field of the helmet dome to form the bright loop of the typical 
loop-like CMEs, followed by the mass of the filament that forms the bright core as 
described in Section 3. When this flux-rope propagates outward, it resembles the 
magnetic clouds that are sometimes observed in interplanetary space as shown in 
Figure 10. We have not been able to determine the threshold energy in this limited 
study because a large number of parametric studies, together with observations, 
are required. This procedure should be considered when an appropriate data set is 
available in the future. 

In summary, a three-dimensional, axisymmetric, time-dependent, magneto- 
hydrodynamic model is, for the first time, presented for a self-consistent description 
of flux-rope emergence from the sub-photosphere. The flux-rope, depending on the 
amplitude of its magnetic energy content, may either remain in place or be ejected, 
either immediately or within the time scale of a day. In both of the two latter 
cases the solar mass ejection consists of two parts: the typically-observed loop-like 
CME and the mass within the flux rope which could be a bright feature following 
the loop as shown in this calculation. However, the present results do not possess 
all the observed features of some global coronal streamer structures that consist 
of prominence, cavity and helmet dome as pointed out by Low and Hundhausen 
(1995); the cavity and dense, low temperature prominence are lacking in the present 
model. We plan to examine this deficiency in die future. On the other hand, we also 
realize that there are still other physical mechanisms such as photopsheric shearing 
(Linker and Mikic, 1995) and heating (Wang et al., 1995) that could destabilize the 
streamer which is worth investigation. To further examine these different scenarios, 
we should look for observational evidences. 
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